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Abstract 

^} \ An exact method for direct calculation of the Jost functions and Jost solu- 

tions for non-central potentials which couple partial waves of different angular 
momenta is presented. A combination of the variable-constant method with 
the complex coordinate rotation is used to replace the matrix Schrodinger 
equation by an equivalent system of linear first-order differential equations. 
Solving these equations numerically, the Jost functions can be obtained to 
any desired accuracy for all complex momenta of physical interest, including 
the spectral points corresponding to bound and resonant states. The effec- 
^) • tiveness of the method is demonstrated using the Reid soft-core and Moscow 

, nucleon-nucleon potentials which involve tensor forces. 
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I. INTRODUCTION 



O I Almost any textbook on the scattering theory has a chapter devoted to the Jost function, 

^ I but none of them gives a practical recipe for its calculation; they provide instead equations 
expressing the Jost function in terms of the wave function. However, to make use of them 
one must find the wave function first which means that the problem is practically solved 
and nothing more is needed. Thus one usually gets a feeling that the Jost function is a pure 
mathematical object, elegant and useful in formal theory, but impractical in computations. 
This is even more pronounced in problems with non-central potentials which couple partial 
waves of different angular momenta. To the best of our knowledge, the classical book by 
Newton is the only one where the coupled partial wave Jost function, which in fact is a 
matrix-function, is considered and the only calculation of the Jost matrix was done in Ref. 
0. The need of such calculations is of course indisputable since many potentials describing 
interactions between molecules, atoms, and nuclei are non-central. 

It is a widely believed that the Jost function for partial waves coupled by potentials 
which are non-zero at the origin, is singular and therefore impractical. However, it is not 
the Jost function but the Jost solution that diverges at r = 0. Of course if one wants to 
obtain the Jost function via the Jost solution, in analogy to un-coupled waves, then the 
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problem of singularity is inevitable and thus one must resort to alternative methods. 

In this work, we show that the Jost function can be calculated directly by simply solving 
certain first order coupled differential equations. These equations are based on the variable- 
phase approach and their solution at any fixed value of the radial variable r, provides 
the Jost function and its complex conjugate counterpart, which correspond to the potential 
truncated at the point r. Such equations were proposed in the past for bound and scattering 
state calculations, that is for calculations in the upper half of the complex momentum plane 

i- 

Exploiting the idea of the complex rotation of the radial coordinate P], r — rexp(i^), 
the method was extended to the unphysical sheet so that the resonance state region was 
also included . This was possible because the coordinate rotation with a sufficiently large 
^ > 0, makes the resonance state wave function quadratically integrable while the energies 
and widths of the bound and resonance states are not affected since the Jost function and 
the positions of its zeros do not depend on r. The effectiveness of the method was demon- 
strated in Ref. for potentials with and without a Coulomb tail where high accuracy for 
very narrow and very broad resonances was obtained. Its suitability in locating Regge poles 
in the complex angular momentum plane was also demonstrated. In this paper we generalise 
the method to non-central potentials that couple different partial waves and show that the 
complex rotation enables us to calculate the corresponding Jost matrix at all points of the 
complex momentum plane of physical interest. In a numerical example we demonstrate that 
a high accuracy can also be achieved. 

The paper is organised as follows. In Sec. II and III our formalism is presented, while in 
Sec. IV the method is applied to several examples and the results obtained are discussed. 
Our conclusions are given in Sec. V. Finally some mathematical details and proofs are given 
in the Appendix. 

II. FORMALISM 

Consider the system of two particles interacting via a non-central potential V{r). Such 
potentials appear in many physical problems as for example in collisions of elementary and 
composite particles with non-spherical molecules and in problems involving spin-dependent 
forces. 

Since this type of potentials are not rotationally invariant, the angular momentum £ 
associated with the interparticle coordinate r, is not conserved. Therefore a partial-wave 
decomposition of the Schrodinger equation results in a system of coupled equations for states 
with different i. In general this system consists of an infinite number of equations, and one 
has to truncate it in order to make it tractable. There are, however, certain problems where 
only few partial waves are coupled to each other, namely, those in which the non-central 
part of the potential stems from the nonzero spins of the particles involved. In such systems 
the total angular momentum, 

J = 1+ s. 
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is conserved and the total spin s defines the maximal number of coupled partial waves by 
the triangle condition 

\J - s\ <i <\J + s\. (1) 

In what follows we will consider a problem of this kind though all formulae remain applicable 
to the more general case of coupled partial waves. 



A. Partial waves for discrete spectrum 



Wave functions "^kjMir) describing bound and Siegert (resonance) states are specified by 
a definite value of the momentum k, the total angular momentum J and its third component 
M, and the parity vr which is omitted in our notation. Such wave functions can be expanded 
in terms of the spin-angular functions 

as 

^kjM{r-) = -Y.y[e]'i^)4](k,r), (2) 

where uj^](fc,r) are unknown radial functions, and [i] stands for the pair of subscripts Is. 
This notation has the advantage that all formulae remain the same for non-central interac- 
tions of spinless particles in which case the symbol [i] stands for another pair of quantum 
numbers, namely, im instead of is. Furthermore all formulae can be generalised to coupled 
hyperradial equations simply by replacing [i] by [L] where L is the grand orbital quantum 
number (for hyperradial equations see Ref. In what follows we drop, for simplicity, the 
quantum numbers JM where this does not cause a misunderstanding. 



Substituting the expansion (||) into the Schrodinger equation, one arrives at the following 
system of coupled equations 



(3) 



where the index [£] runs over all combinations of £ and s allowed by the triangle condition 
(0) and the parity conservation law. The elements of the matrix W are those of the operator 
V{f) (we use ^=1), 



tJM\ 



(4) 



sandwiched between the spin-angular functions (m is the reduced mass). We assume that 
these matrix elements are less singular at the origin than the centrifugal term. 



limr2vr[^][,,](r) 



0, 



and vanish at large distances faster than the Coulomb potential, 

0. 



lim rWro^,m[r] 



(5) 
(6) 
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B. Partial waves for continuum 



The scattering state wave function "^j^g^ir) is defined by the real vector k (the momentum 
of the incoming wave), total spin s, and its third component fi. The partial wave decompo- 
sition for this function is more complicated than for the bound and Siegert states because 
the scattering state depends on the direction of the incident momentum k. Since "^j^s/ii^ 

depends on the two vectors k and f, we have to perform partial wave analysis in both the 
momentum and coordinate space 



* " jM[e']e 

where 

^Mf(k)=^-'EC/i>^.™(k), 

m 

and the radial wave function r) obeys the same equation as Mj£](/c,r), 

li"] 

Physical solutions of Eqs. (0) and (|^) are defined by the requirement that they must be 
regular at the origin, 



(9) 



and by certain, physically motivated, boundary conditions at infinity, 

(10) 

which are different for the various problems (bound, scattering, and resonant states) and 
change drastically when we go over from one case to another. It would be, of course, more 
convenient to deal with a universal boundary conditions imposed at a single point. This 
can be achieved if we consider the general regular solution of Eq. (^ defined only by the 
condition (H) and not subjected to any restrictions at large r. This is described next. 



C. Regular basis 

From the vast variety of solutions obeying (P) and having different behaviour at large 
distances we choose only those which are linearly independent. They form the fundamental 
system of regular solutions which we call regular basis. Any regular solution with a specific 
behaviour at large distances, is simply a linear combination of the basic solutions. Thus, 
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instead of having different mathematical procedures for the various types of physical prob- 
lems we can have only one for the regular basis. In the next section we show how the bound, 
scattering, and resonant state wave functions can be constructed in terms of such a basis. 

Let us consider Eq. as a matrix equation. Then, each of its solutions is a column- 
matrix. From the general theory of differential equations is known that there are as many 
independent regular column-solutions of Eq. as the column dimension, i.e., the number 
of equations in the system. These columns can be combined in a square matrix 
satisfying 

dl + e- i{i + l)/r^] <^m'ik,r) = Y,W[eW']{r)<^[e"]m{k,r) , (11) 



with 



^im']ik,r)^^0, V (12) 



Since Eq. (^TJ) is of second order and singular at r = 0, the condition ([T^) cannot be reduced 
to the simple requirement 

<l>[,][,,](A;,0)=0, V 

as the behaviour of each element of the matrix ||$(/c,r)|| and its first derivative (9r.||$(A;, r)|| 
in the immediate vicinity of the point r = are also needed. 

Like in any other space, the basis can be chosen in an infinite number of ways by speci- 
fying the behaviour (|12D. The possible choice of the condition ([l2|) , however, is not entirely 
arbitrary. It was shown in Ref. [jlO[ (see also an alternative proof in the Appendix) that for 



potentials fulfilling the restriction (|^), the regular columns are linearly independent only if 
they vanish near the point r = in such a way that 

$M[£'](fc,r) 

S — — -^m']- (13) 

The primed angular momentum in the denominator means that in each row of the matrix 
the elements, situated further to the right, vanish faster when r — > 0. Even with- 
out a formal proof it is clear that such condition guarantees the linear independence of the 
columns since in each row the elements have different behaviour (different power of r) near 
r = 0. If instead of the rows, we look at the columns, we get from Eq. (|T3p that in each of 
them all off-diagonal elements are infinitesimal as compared to the diagonal one. 



Yet in spite of the restriction ([T3|), we still have some freedom in specifying the deriva- 
tives dr^[e][ii]. Indeed, we can choose at least an arbitrary r-independent coefficient in each 
element of the matrix || ||. We mention here that it was this freedom which was ex- 
ploited by Newton in his procedure of regularisation of the integral equation for [0]. 

To be consistent with the definition of the regular solution for the un-coupled partial 
waves |TTI], we choose the normalisation coefficients in such a way that 
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where je{z) is the Riccati-Bessel function |jT2[. This condition defines the leading terms of 
the near-origin behaviour of the diagonal elements of the matrix || $[£][£/] || and their first 
derivatives. The off-diagonal elements, however, remain obscure since the condition (|14D 
only implies that ~ o{ji'),i ^ £'. As pointed out by Newton [Q], it is impossible to 

define them unambiguously by boundary conditions which are independent of the behaviour 
of the potential near the origin. In the Appendix we show that this is indeed the case and 
we give a simple recipe for obtaining series expansions of all matrix elements of r)|| at 
r ~ 0. Of course, the terms of such series depend also on the potential. 

In principle, from the knowledge of such expansions, we could calculate the matrix 
||$(/c,r)|| by solving the Schrodinger equation ( ]TT| ) directly. It is, however, much more 
convenient to transform Eq. (|lT]) into another equivalent form suitable for a construction of 
different physical solutions. For this we introduce two new unknown matrices r)|| 
and assume the following ansatz for the regular solution 

^[emik,r) = i (A:r)4+J,,(/c,r) + h\--\kr)F{~^,^ik,r)] (15) 

where the Riccati-Hankel functions h^^^ are linear combinations of the Riccati-Bessel and 
Riccati-Neumann functions h^^\z) = je{z) ± ini{z) |12[. The reason for choosing such an 
ansatz will become clear in the next section. Here it is sufficient to say that the explicit im- 
plantation of the functions h^^\kr) into the construction of the basis guarantees the correct 
asymptotic behaviour of the basic solutions at large r. 

Since instead of one unknown matrix ||$|| we introduced two matrices ||-F''^^||, they must 
be subjected to an additional constrain. The most convenient is the Lagrange condition 

/il+^(A:r)a,4+],](fc, r) + /i^)(fcr)a,4nVfc, r) = , (16) 

which is a standard choice in the variable-constant method for solving differential equations 
p!3| . Substituting (p!5D into Eq. (|Il|) and using the condition (p!6|), we obtain the following 
first order coupled differential matrix equations 



drF[^^,^{k,r) = ± [h\t\kr)Fl;,l^ik,r) + h\-\kr) F^-^^^.^ik , r)} . 



,(T), 

2ik 

(17) 

As we show in the Appendix, the boundary condition (p^ can be rewritten as 



lim 



je{kr)Ff^^^^,^{k,r) 
ji'{kr) 



S[i][e'] , (18) 



and a series expansion for ||Fj^j'|],j(/c, r)|| near the point r = can be found iteratively as 
follows 
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Fim\k^r)=6[e][n, (19) 



(20) 



/i^/, '*(A;r)F[^jj]^j ■'(/c, r) + /i^„-*(A;r)F[^jj],*j \k,r] 



dr . 



where the arbitrary constant for the indefinite integral is assumed to be zero. We shall give 
an explicit example for this expansion in Section IV. 

In contrast to Eq. (|TI|) the equations for the new unknown matrices ||F'^^)(/c, r)|| are of 
first order. However the point r = is generally a singular point because the Riccati-Hankel 
functions have the short range behaviour ~ and according to (|^) the potential may also 
behave near this point as ~ r~^^~^\ e > 0. Therefore Eqs. (p!7| ) cannot be solved with the 
boundary conditions for ||F*^^)(A;, r)|| at r = 0. Instead, we may use the analytical solutions 
of them in a small interval (0, 6] (in the form of the above series expansions) and impose the 
boundary conditions at r = 5. 

Certain elements of the matrices ||F*^^^(A;, r) || could diverge when r ^ (see the Ap- 
pendix). This, however, does not cause any problem in the iterative procedure (^) since the 
integration is represented by an indefinite integral. It is emphasised that although certain 
elements of HF*^^)!! diverge, the matrix ||$|| always remains regular because at small r the 
matrices converge to each other, that is they are transformed into some matrix ||y4|| 

which is the same for both of them, H-F^^-* || \\A\\ (see the Appendix), and the singularities, 

if any, are compensated by the behaviour of the Riccati-Bessel function je{kr), 

[^^^^(^0 + hi-\kr)] A[,-^^e]{k,r) = j,{kr)A[,^[e]{k,r) . (21) 

The system of equations ( p^ together with the boundary values F[^^],](fc, S) ~ Fl^^^i^\k, S) 
represent a well-defined differential problem, a solution of which gives the regular basis in 
the form (|TB|). 

D. Jost matrices 

It can be proved (see the Appendix) that for Imk = the right hand side of Eq. (|l^ van- 
ishes when r — i> oo. Since the corresponding derivatives also vanish, the functions Fj^^],j(/c, r) 
become r-independent and thus, for momenta corresponding to the scattering states, we have 

(k, r) \ (fcr) (A:) + h^f\kr)T\^Uk)\ (22) 



where 



^m']^k) = \m^F[^,^{k,r) . (23) 
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We may call these r-independent matrices ||^*-^''|| as Jost matrices and the products 
h^[^\kr)F^^^,^{k^r), which behave asymptotically like ~ e^*^'", as Jost solutions. 

For complex values of k the above limits generally exist in different domains of the 
complex /c-plane, namely, in the lower half (Imfc < 0) while in the upper 

half (Im/c > 0). This is because, according to Eq. (p!7|), the derivatives drF^^^k,r) are 
proportional to hY\kr) with 

hf\kr) ^-^^iexp [±i{kr - ^^^/2)] (24) 

vanishing in different domains of the fc-plane, namely, 

ht\kr)^Q, Im(A;r)>0, (25) 
h[~\kr)^Q, Im(A;r)<0. (26) 

Thus, in general, the only area where the limits (|23| ) simultaneously exist is the real axis Q. 
However, for a particular class of short-range potentials (decaying exponentially or faster) 
the upper bound for the existence of JF(+) is shifted upwards and the lower bound for J^^~^ 
downwards, which widens their common area to a band. 

The difficulty concerning the existence of the limits (|23|), can be circumvented in the 
same way as for central potentials |§|J^. Indeed, the conditions ( p5D and ( p6D involve the 
imaginary part of the product kr and not of the momentum alone. Therefore, if, for example, 
Im {kr) is negative we can make it positive by using the complex rotation method which we 
describe next. 



E. Complex rotation 

In this method the radius r is replaced by a complex one, viz., 

r = xexp(i^), a;>0, < |0| < ^ . (27) 

The idea of complex rotation of the coordinate is not new. Many years ago, during World 
War II, Hartree and his co-workers at Manchester university used such rotation to solve 
certain differential equations describing radio wave propagation in the atmosphere (for more 
details see Ref. |]14|)- Nowadays the complex rotation is widely used for locating quantum 



resonances by variational methods 0. In contrast, our equations together with the complex 
rotation can be used to locate resonances in an exact way. 

Applying the complex transformation (^) to Eq. (p!7|), we obtain 



It can also be proved that both these hmits exist at all spectral points corresponding to bound 
and resonance states (see the Appendix) 
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ri (28) 



The purpose of the rotation (|27|) is to make the imaginary part of the product kr positive or 
negative, in calculating JF(~) or JF(+) respectively, at points on the fc-plane we are interested 
in. Thus we have Im fcr > for all points above the dividing line shown on Fig. 0. This line 
defined by the negative angle 6 in the /c-plane results from the rotation ( ^Tf ) in the r-plane 
with positive 6. 

If the potential matrix ||W^(r)|| is an analytic function of the complex variable r and 
obeys the conditions |^) along the ray (0), then the limit 

il^4i[i'](^'^^^')=-^[S'l(^) (29) 

exists and is finite for all k on and above the dividing line [— cxoe"*^, +ooe~*^] (for the relevant 
proof of this statement see the Appendix). At the same time the limit 

exists and is finite for all k on and below the dividing line. Moreover, when the limits (^91) 
and (|30|) exist the values of ^[^^ii-^{k) are independent of the rotation angle 6 as the Jost 
function is i — independent and hence 6'-independent. Thus, the limits ( p9D and ( pO]) give 
a unique analytic continuation of the Jost matrices to the lower and upper halves of the 
complex fc-plane respectively. To calculate the ^[£]j]/](^) for Imfc < we need to solve Eqs. 

(p8|) at a sufficiently large positive 6, and the ^^£^^p^{k) for ImA; > at a sufficiently large 
negative 6. 

We note that though the ansatz (^) is suitable for large distances (see forth the next 
section), it is not good for numerical calculations in the vicinity of r = 0. Indeed, near this 
point the singularities of h^^\kr) and h[ \kr) are cancelled. Although this does not cause 
formally any problem, in numerical calculations the cancellation of singularities is always 
a source of possible numerical errors. These errors increase with increasing ^ since in this 
case becomes more singular. Therefore, for larger £ the point r = 5 must be shifted 



further from the origin. This in turn, requires more iterations of Eq. (|20D to obtain the 
boundary values F^^\k,6) to a required accuracy. 

Eq. (|2T| ) hints to another way to handle numerically the boundary condition problem. 
Since ^ ^ 

- {hj + hj)= and —{hi - hj) = , 

we may introduce a new pair of matrices 



Aye^[e^{k,r) = i [F[^J|],](/c, r) + FjijjJ,j(fc, r) 



2 



B[i][i']{k,r) = — \Fr[tl{k,r) - FLl{k,r] 



(31) 



2i 
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which transform the ansatz ([T5|) into the form 

^[i][e']{k,r) = ji{kr)Aii][e>]{k,r) - ni{kr)B[i][(.>]{k,r) , (32) 
and the corresponding hnear combination of equations (0) gives the alternative form 

ni{kr) 



drA[(.][p]{k,r) 



"1 r 



[£"] 



je'{kr)A[e"][e']{k,r) - ne"{kr)B[g»][i>]{k,r) 



jt{kr) 



Yy^[m"]i'^) h"{kr)A[ej,^^[e]{k,r) - ne"{kr)Bii»][i>]{k,r) 



Likewise, the iterative procedure ([T^, [20| ) transforms into 



(33) 



(34) 



X 



>/(A;r)A[^"''j[^,](A;,r) - n^//(A;r)5[^"/j[^,j(A;, r)J c?r , 



in) 



(n) 



dr . 



(35) 



The representation of ||$|| in terms of and is equivalent. From a practical 

point of view, however, it is more convenient to start the integration of Eqs. ( |33D from the 
boundary values A[£][£'](fc, (5), (/c, 5) and at some intermediate point Tint (far enough 

from the origin) to go over to Eqs. (|T7D with starting values Fj^j^],j(A;, rj„t) obtained from 
y4[£][^/](A;, rj„t) and -B[£][f'](fc, Tj^i) via the linear combinations (|3T1) . 



One may argue that we can abandon equations ( [17|) altogether and integrate instead Eq. 
(0) on the whole interval [5, rmax]- However, ||F^"''''(A;, r) || and ||F^~-'(/c, r)|| have finite limits 
(r oo) in different domains of the complex A;-plane (below and above the dividing line 
respectively). The only points where they have limits simultaneously are the spectral points 
and the dividing line itself. And since \\A\\ and involve both they have limits 

only at these points. Therefore, to obtain the Jost matrix we should start at small r with 
Eq. (|33[) and finish at large Vmax with Eq. (p!7|). 



III. PHYSICAL SOLUTIONS 

In what follows we shall describe how we can obtain a physical solution from the regular 
basis. In general, each column representing a physical solution is a linear combination of 
the basic columns. 
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[£'] 

(36) 

with the coefficients ||c|| defined by the physical boundary condition ([l0| ) at large distances. 

A. Bound states 

The bound state wave function vanishes at large distances as 

5] $[,][^,] (A;, r)c[,,] %e-l'=l^ , 

where N^e] are the asymptotic normalisation constants. In this equation the function 
(A;, r) can be replaced by its asymptotic form (p2D, i.e., 

1 E K^(fcr) J-g,j(A:) + /^J-)(A:r)^^(fc)] q,,] ^^^O . (37) 
[€'] 

For bound states Imk > and the Riccati-Hankel function h^^\kr) decays exponentially 
while h[~\kr) grows exponentially. Therefore the condition (^) can be fulfilled only if we 
find coefficients C[£] such that the diverging functions h[~\kr) of different columns cancel 
out, that is, if 

E-^S'](^)c[^']=0- (38) 

This system of homogeneous linear equations has a nontrivial solution if and only if 

det||J^(~)(A;)|| =0. (39) 

Therefore, we can locate all possible bound states by looking for zeros k^ of the Jost-matrix 
determinant on the positive imaginary axis (see Fig. 0). For each zero ko thus found, 
the coefficients cj^j are then uniquely determined by the system ( P8| ) apart from a general 
normalisation factor which is finally fixed when the physical wave function, 

^fcoiM(r) = 1 E y^'i^) [hi'\kor)Fli},^iko,r) + h[-\kor)Fl^^,^iko,r)] q,] , (40) 
me'] 

is normalised. The contribution of each element of the column 

Mko, r) = W [h[-'\kor)Flg,^{ko, r) + h^f \kor)F{-^,^{ko, r)] q,,] (41) 

into the normalisation integral represents what is usually called the percentage of the cor- 
responding partial wave. 
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B. Scattering states 



The scattering states normalised to the 5-function, 



are defined by the following asymptotic condition 



{2n 



1 3/2 



ikr 



e'^^'^^Xs/. + J2 is',,'s,.{k-, k)x. 



S'fJ.' 



(42) 



where is'n'sfiik', k) is the scattering amplitude. The partial wave decomposition of Eq. 
gives for the boundary condition ([T0|) (see Ref. [IH) 



umm{k,r) 



1 r 



r— >oo 2 



h\,\kr)5^ni] + hV\kr)S{,,^^,^{k) 



(+)< 



(43) 



Therefore the choice of our ansatz ([T5|) for the regular basis is natural and suitable not only 
for constructing the bound states but the scattering states as well. Indeed, comparing 
with (p2|) we find that the coefficients C[^/][^] in (|3^) should be chosen as 



C[e'][e] = \\^ie^]li]{k)\\ ^ , 
which gives us the S-matrix in the form 

||5(A:)|| = ||^(+)(A:)||.||^(-)(fc)||-^ 



(44) 



Thus, the normalised scattering wave function can be constructed from the regular basis as 
follows 



2TTkr 



EEE 



(45) 



JMe [£'] [£"] 

3^[''f {h\t\kr)Fl;l,^ik,r) + /^(r)(A:r)FH,„,(A:, r)} WJ^^^.^m-' ■ 

The scattering phase shifts together with the mixing parameters can be found from the 
S-matrix given by Eq. 



C. Resonances 

The resonance (or Siegert) states are described by wave functions which at large distances 
have only outgoing spherical waves 

^UMir-) = I E 3^[r(r)nf.](A:, r);:^-"^. (46) 

More precisely, the resonant boundary condition at large r is 
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ui^{k,r) ^-^Nf,^{k)hi-'\kr), (47) 

where N^^-^{k) are the partial wave normalisation constants. Since there is not a generally 
accepted convention about the normalisation of Siegert states, the N^^-^{k) may involve an 
arbitrary coefficient. Solutions of the Schrodinger equation, with the long-range behaviour 
(^), may exist only at discrete points of the complex /c-plane, situated below the real axis 
(see Fig. The corresponding radial wave functions are regular at the origin and are 
therefore linear combinations of the regular basis 

r) = lj: {hi'\kr)Flg,^ (k, r) + h'~r\kr)Fl^^,^{k, r)] cy,^ . (48) 

This equation can be used to locate resonances. Similarly to the bound states, we simply 
require that 

m 

which means that U[e]{k,r) is expressed only in terms of h^i^\k,r) at large distances. Since 
6 can be chosen such that Im kr > 0, the Riccati-Hankel function h^^\kr) decays exponen- 
tially when \r \ oo. This in turn means that the Siegert states are quadratically integrable 
according to (|47|) and hence they acquire the same properties as the bound states. There- 
fore, with the coordinate rotation the bound and resonance states can be treated in the same 
way. In particular, the position of a resonance is defined by Eq. (p9D, and its (rotated) wave 
function by Eq. (|40| ) where r is now complex. 

Since this wave function is square integrable we can normalise it to unity and thus the 
normalisation constants N^^-^ in (^) can be fixed in a natural way. Therefore the proposed 
method enables us not only to locate the position of the resonances as zeros of the Jost 
matrix determinant but also to obtain the correct normalisation constants N^^^^. In order 
to obtain physical (unrotated) Siegert wave functions, the Schrodinger equation must be 
integrated along real r inwards using the boundary condition (^), the found momentum /cq 
and normalisation constants. This integration will provide automatically a wave function 
which is zero at r = because ko is a spectral point. 



IV. EXAMPLES 



In order to demonstrate the effectiveness of the method we consider as as example the 
nucleon-nucleon (NN) interaction in the triplet spin-state, i.e., when the total spin s = 1. 
From the Pauli principle it follows that s is conserved (see for example Ref. |]15|); and thus 



the sum over [i] is reduced to Y.e- The triplet NN potential can couple only at most two 
partial waves, with i = J —1 and i = J + 1, a.s the state with i = J has different parity and 
therefore must be excluded. We consider here the even state of two nucleons with J = 1, in 
which they can form the deuteron. The partial wave decomposition of this state consists of 
coupled S and D waves. 
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The corresponding NN potential includes the following three most important terms 

V{f) = K(r) + Vt{r)S^2 + V,.s{r){I- s) , (50) 

known as the central, tensor, and spin-orbit potentials. The second term contains the tensor 
operator 

3 

•^12 = • ^)('?2 ■ r) - (cti ■ 0^2) 

which is responsible for the coupling of different partial waves. This is clearly seen from the 
structure of the matrix \\W\\ defined by Eq. (H), which in the case of s = 1, J = 1, and 
TT = +1 reads (see for example Ref. |l|]) 



/ Woo W02 \ 




( 


2V2Vt \ 




= 2m 






V W20 W22 ) 




[2V2Vt 


K - 2Vt - 3Ves ) 



Wu'{r) — > au'T ^ + hu, + c«/r + . . . 
1 — >o 



(51) 

where the subscripts of correspond to the two partial waves £ = and £ = 2. 

In accordance with the restriction (^ we may use at short distances the series expansion 

(52) 

which is general enough to include potentials having a soft core (a«/ 7^ 0). 

Using ( |52D we can obtain the corresponding expansions for Au/{k, r) and Bur{k, r) via the 
iterative procedure (|3^ , |35|) which are needed to impose the boundary conditions Aii/{k,S) 
and Bui{k,6) at a small 6. These expansions near the point r = can be obtained in an 
explicit form. For this we replace the Riccati functions in the indefinite integrals (|35D by 
their series expansions 

j^ikr) = kr - + ^ ^ 



j2{kr) 
no{kr) 
n2{kr) 



6 120 
(kr)^ (kr)^ (krY 
"TS 210" ^ 7560 
{krf {kry 



-1 . 

2 24 

3 1 {kry 

' {krf ~2 8~ 



+ ■ 



Substituting (|34D into (^) and performing the integrations over r for each element of the 
matrices 11 A 11 and we obtain 



||A«(fc,r)|| 



1 + floor —r^r^ * 
45 

3a20 -1 T , ^22 



|i?«(A;, 



/ aook 



020 A;^ 
60 ' 



Q-02^'^ 4 \ 



60 

022 fc^ 



1350 / 
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where we retain only one additional term in each matrix element. For the second iteration 
we get 

/ 1 I I Coo 2 

1 + aoor + — r 



\\A('\k,r)\\ 



ao2^^ 3 , 
-r + 



45 



60 



3^20 -1 , 3^20 , 

r H — 77;- mr 



V A;2 
/ ooo^ 

020^^ 



k^ 

kCoo 3 
— 



1 022 , C22 2 

1 H r H r 

5 10 



V 



60 



20 5 
r 



where the constants (u' are defined as 



/ 




^00 , t 

^ + ^00 



a2oCtoo 



V 



'20 



75 

Q02Q20 
4 

_ Ct22a20 



ao2k''^ 
60 

a22k^ 
' 1350 ' 



A;'Co2 5 \ 

r 

75 



22 7 
r 



1575 / 



'^00^^02 , 7 , Cf02fl22 \ 
+ O02 + 



12 

Ct20Ct02 
12 



6 



fe22 + ^ 



6 / 



The iterations can be continued in the same manner with each iteration adding a new term 
to each matrix element having a higher power of r than the previous one. All matrix ele- 
ments of ||y4|| and ||i?|| are regular when r —>■ except for the left bottom corner element of 
\\A\\. It has two singular terms, ~ and ~ Inr. The next iteration however gives for it a 
vanishing term of the kind ~ rlnr, and all subsequent terms are vanishing even faster. 

The above expansions illustrate the fact that at small r the matrices and ||-F*^~)|| 

converge to each other. Indeed, since \\B\\ is infinitesimal as compared to \\A\\ the second 
term in the linear combinations = ||y4|| ± i\\B\\ vanishes when r 0. 



The singularities of the matrix r)|| refiects the main difficulty of the theory 

of coupled partial waves, which precluded its development in the past. As we show in the 
Appendix, the matrix elements of r)|| with i > i' are always singular at r = if 

the potential is nonzero at that point. Despite this the functions r) remain always 

regular because in Eq. (|32D the singularity of r) is compensated by ji{kr). In our 

case we have 





(53) 



and therefore the singular term A2o{k,r) appears only in the product j2(^'")^2o(^? which 
vanishes at the origin as ~ r^. Substituting the series expansions of all functions involved 
into Eq. (ESI) we obtain 



15 



kr H r + 



2 I CoO^ 3 , 4i \ 

I r +C(r In rj 



6 



180 



^4 + ^ (^(^6) 



300 



(^20^ 2 I C20^ 



r Inr 



fc3 
15' 



(54) 



\ 4 ■ ' 5 ■ 25 ■ ' 15' ' 90 ■ ' 210 

which exphcitly demonstrates that our regular basis obeys the boundary condition (|I^) and, 
as mentioned in Sec. II, in each row the convergence to zero increases from left to right while 
in each column the diagonal elements have the lowest vanishing speed. The linear indepen- 
dence of these columns is also apparent. As far as the off-diagonal elements are concern, even 
the leading terms of them depend on the behaviour of the potential, i.e., it is impossible to 
specify the boundary condition for them in a general form independently from the potential. 



Using a quite different approach, Palumbo in Ref. derived recurrence formulae for 
direct constructing the series expansion of the regular basis for the potentials of the type 
(^). Having performed few iterations of the Palumbo's formulae, we found that they gen- 
erate the same terms which are given in Eq. (0). This is yet another confirmation that our 

I) and flSl, BSl) are correct. 



iterative procedures (|l£ 



In order to test the proposed method numerically, we chose two different NN-potentials 
of the type (|52D. The first one is the Reid soft core potential (RSC) [0, which for J = 1, 
s = 1, TT = +1 has the following form 



-2ar 



-Aar 



-Gar 



K(r) 



Vt(r) 



ho- 



ar 



+ hi- 



+ h2 



ho 



+h4 



+ 



ar 



-4ar 



ar 

3^3 
[ary {ary 



ar 



+ h3 



ar 
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(ar 



+ 



[ar] 



ar 

-4ar 



+ h. 



ar 

-6ar 



VfJr) 



ar 



+ h-j 



ar 



with 



ho 
h3 
he 



-10.463 MeV 
9924.3 MeV , 
708.91 MeV, 

The second potential used, is the Moscow potential |jl7|, 



hi = 105.468 MeV, 
/i4 = 351.77 MeV, 
hr = -2713.1 MeV, 



/i2 = -3187.8 MeV ^ 
/i5 = -1673.5 MeV; 
a = 0.7 fm"^ 



K(r) = Vie-^"-" +¥2(1- 



--yr 



(3r 



Vt{r) = V2 



(3r (/5r)2 



'■yr 



3 e" 



-l3r 



pr 



(55) 



(56) 
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with 

= -466.74 MeV , V2 = -10.69 MeV , 
/5 = 0.6995 fm-\ 7 = 3fm-\ r] = 1.6im~^ . 

The RSC potential has a strong repulsion at small distances, 

2m / ho + hi + h2 + /i3 2v/2(/i4 + /is + 23.5/io) \ 

V V 2v^(/i4 + /i5 + 23.5/io) /ii + /i2 + /is - 2/i4 - 2/i5 - 46/io - 3/i6 - 3hj ) 

_ 2m / 6832 MeV -4434 MeV \ 
~ ~ 1^ -4434 MeV 15979 MeV ) ' 

In contrast the Moscow potential has very strong attraction instead and sustains, apart from 
the deuteron bound state, a very deep bound state known as Pauli Forbidden State (PFS). 
In this case the expansion (|5^) begins from the second term, 

-513 MeV -7156 MeV 



2m 



2m 



-7156 MeV 4547 MeV 



Both potentials describe the deuteron properties and the np-scattering quite well despite 
their completely different short range behaviour. 

To begin with we consider real energies corresponding to bound and scattering states. 
We integrated Eqs. (0) by the Runge-Kutta method from r^m = 10~^ fm to Tint = 1 fm 
with the boundary conditions HA*-^-* (/c, rmm) || and ||i?^^)(/;;, rmm)||- Then from r^nt = Ifm we 
integrated Eqs. (0) up to r^ax = 20 fm where the functions F*^^^(/c, r) attain their limits 
(^). Repeating such calculations with different values of the momentum k corresponding 
to points on the positive imaginary axis, we found that the equation 

det||F(-)(A;,w)|| =0 (57) 

is fulfilled at the points ko given in Table I. The binding energies and the percentages of the 
D-waves {D%) are also given in this table. For comparison Table I contains the energies 
and D% obtained originally in Refs. |T^-[T^ by the authors that constructed these poten- 
tials. The other observables such as the mean square radius and the electric quadrupole 



moment of the deuteron are also the same as given in Refs. |lT6| , |r7| . Due to the presence of 
the deep unphysical PFS state in the Moscow potential the deuteron state is an excited one 
and therefore its wave function has a node at Tc ~ 0.59 fm. 



For real positive k, i.e., for scattering states, we performed calculations with the same 
rmin, rint, and Vmax- According to Eq. (||) the product 

\\F^^\k,r^,ax)\\ ■ \\F^-\k,rmax)V = \\Sm 

gives us the S'-matrix which contains information about the scattering observables. The 
nucleon-nucleon S'-matrix usually is parametrised in terms of the so-called bar phase shifts 
and mixing parameter, introduced in Ref. [0, as follows 



[ "S'oo 


S'02 




\ 'S'20 


S'22 J 





e^'~^° cos 2e ie'^^o+h) 2e 
^gi(5o+S2)gjn2e e2*^2cos2e 
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^From this matrix equation one gets 



So 

52 



\nSoo 


/- 


C2 
^02 


2i \ 


'S'oO'S'22 


In 5*22 


/- 


C2 
^02 


2i \ 


5'oO'S'22 



e = — arccos 



'S'oO'S'22 + 



4 \5'oo>S'22 — >S'o2 / 

The obtained phase shifts and mixing parameters are in agreement with the values given 



in Ref. ||T6[, for all collision energies examined (up to Ecm. = 176 MeV). For example, at 



Ecm. = 12 MeV {k = 0.53793 fm"^) the RSC potential givesQ 
||FW(fc,r 



11 _ / 0.16247 X 10" T^O. 18386 X 10^2 -0.60905 x 10^ ± zO.68923 x 10« 
, ' max) II - I _o. 13390 X 10^3 ± i0.61083 x 10" 0.50194 x 10^ T «0.22898 x 10^ 



and 

II c(i.) II _ ( -0.95640 + i0.28507 -0.06232 + i0.01229 \ 
" ^" ~ 1^ -0.06232 + i0.01229 0.99303 - z0.09933 ) 

^From the latter 5'-matrix we obtain the following scattering parameters 
So = 1.4288 , ^2 = -0.04995 , 2e = 0.06357 , 



which practically the same with those obtained by Reid ||T6| 



60 = 1.426 , 52 = -0.050 , 2e = 0.064 , 
via direct a solution of the Schrodinger equation. 



The Moscow potential at the same energy gives different Jost matrices 
\\F^^\k,r^a.)\\ = 



361.16 ±i4347.4 -0.08458 ± il. 0183 
30864^^1423.1 7.2291 ^^0.33331 



but practically the same S-matrix 
11-5(^)11 = 



-0.95557 + iO.28772 -0.06283 + i0.01252 \ 
-0.06283 + Z0.01252 0.99268 -^0.10059 ) 



and therefore the same phase shifts and mixing parameter 

5o = 1.4275, ^2 = -0.05059, 2e = 0.06411. 



^The calculated H-F^^-* (A;, rmax) || and ||F( \k,rmax)\\ are complex conjugate to each other at least 
within five digits 
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The huge values of the above Jost matrix elements are due to the behaviour of the potentials 
at small distances, which generates large values of the derivatives drWF^'^^k, r)\\ that pushes 
up the absolute values of the functions r). To demonstrate an opposite example, we 

solved Eqs. (|1^) at the same energy {Ec.m. = 12 MeV) with a rudimentary nucleon-nucleon 
potential (also in the deuteron channel) consisting of two Yukawa-type terms 
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V{r) = ,/"P(-"^/P^) + ^^exp(-^rM)^^^ 



r/pc 



r/pt 



(58) 



= -22.7 MeV, Vt = -W.9MeV , 
u = 2.l2, pe = 2.47fm, pt = 3.68fm. 

Like the RSC it is singular at r = 0, but the coefficients 



\0'U' 



/ -22.7 MeV -45.9 MeV \ 
^"^/'c _45.9MeV 9.8 MeV ) 



for its expansion (^) are in two orders of magnitude less than those for the RSC potential. 
As a result the Jost matrices calculated at Ecm. = 12 MeV are 

,, 0.27328 ±z0.87845 0.11204 ± iO. 23730 \ 

" ^ '"""^""^'l ~ 1^ 3.8437 T«0.37653 1.6876 ^ i0.15658 J ' 

while the S'-matrix is 

II c(i,\ II _ ( -0.98755 + iO. 10546 0.11564 - i0.015980 A 
" ^" ~ 0.11564-^0.015980 0.97914-^0.16631 ) 

are 5o = 1-53, 82 = -0.0847, and eo = 0.117. 



So far we dealt with momenta on and above the real axis (Im k > 0) and therefore the 
coordinate rotation (|27| ) was not needed. Consider now a point k which is under the real 
axis. To obtain the Jost matrix ||jF(~)(fc)|| in this domain of the fc-plane, we must integrate 
the rotated equations ( P^D (with > 0) along the real variable x = \r\. Similarly to the 
case with ^ = 0, in the immediate vicinity of the point a; = it is convenient, for numerical 
reasons, to replace by their linear combinations (^Tj) and the equations ( |25D by the 

corresponding linear combinations of them. The resulting rotated equations for \\A\\ and 



||-B|| as well as the rotated iterative equations follow immediately from (|33D , (|3^ , and (^ 
after simple replacement of r by xexp{i6). 



To demonstrate the ability of the method to deal with momenta of the fourth quad- 
rant of the fc-plane, we calculated the Jost matrix for the RSC potential at A; = 
0.5 exp(— z0.37r) fm~\ i.e., at Ec.m. ~ (—3.20 — 29.86) MeV, with three different values of the 
rotation angle: ^ = 0, 0.357r, 0.47r. The results obtained are given in Table II. The first line 
of this table demonstrates that the unrotated equations cannot give a correct 6'-independent 
Jost matrix when Imfc < 0. The matrix obtained with = is also r-dependent, i.e., has 
no limit when |r| 00. If however 6 is large enough, such that Imfcr > and the point k is 
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above the dividing line, then \\F^~\k, XmaxG^^)\\ does not depend on 9 (compare the second 
and third hnes of the table) and Xmax- To achieve the Xmax^independence when 6' 7^ 0, we 
have to go further afield because the potential vanishes along the xexp{i6) slower than along 
the real r. Thus, the results given in Table II were obtained with Xmax = 50 fm. 

The number of digits which are unchanged under the rotation show the accuracy achieved. 
In this connection it is interesting to note that the accuracy of the second column of the Jost 
matrix as well as of the determinant turn out to be always much higher than that of the Jost 
matrix first column. This is exemplified by the two last lines of Table II. It is interesting to 
know that the correct value of the determinant can be obtained even with crude boundary 
conditions and large tolerance of the Runge-Kutta procedure. Thus, the correct binding 
energy of the deuteron (given in Table I) can be obtained even with boundary conditions 
11^^^) (A;, r^m) II and \\B^^\k,rmin)\\ when the first column of \\F^~\k,rmax)\\ is even wrong. 
The only explanation for this is that somehow both elements of this column get the same 
erroneous term which cancel out in the determinant. This observation means that the pro- 
cedure of locating spectral points is less demanding and less delicate than the calculations 
of the corresponding wave functions. 

None of the potentials (^,^,^) generates resonances (at least at reasonably low en- 
ergies). To the best of our knowledge, non-central potentials generating resonances have 
not been published yet. Thus, in order to demonstrate the ability of our method to locate 
Siegert states we constructed an artificial potential with a rich spectrum. For this we used 
the well-known central potential 

\4(r) = 7.5r^exp(-r) , (59) 

which is widely used as a testing case for new methods that locate resonances (see for ex- 
ample Refs. |[7|j21|-p3|) . It is usually assumed that this potential is given in atomic units. In 



order to be consistent, however, with the potentials used in this work, we assume that it is 
given here in MeV. Then for h'^/2m = l/2MeVfm^, the numerical values of the resonance 
energies are the same in MeV and in atomic units. 

This potential generates a sequence of S-wave resonances (see Ref. 0) which cannot 
be significantly displaced if we add very weak interaction in the D-wave and a weak S-D 
coupling. Thus the non-central potential coupling of the S and D partial waves which is 
represented by the matrix 

\\Wuir)\\ = 2m ( 7-5r^exp(-r) -Ar^exp(-r) \ 

II «v ;ii \^ _;\^^exp(-r) -Ar'^exp(-r) y ' ^ ^ 

should generate resonances at least when A is small. Bound states can also be generated 
by increasing A since we have chosen a negative sign for the D-wave potential and for the 
off-diagonal elements of the matrix (pO|). 



Of course the use, as a testing case, of a potential with a weak D-wave and weak coupling 
is rather undesirable. Therefore, we gradually increased A up to the value A = 15 MeV. The 
resulting spectrum is given in Table III and depicted in Fig. ^ It consists of eight bound 
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states and a sequence of resonances (we show only six of them, nearest to the real axis). By 
the large open circles in Fig. ^ we also show the positions of the first three resonances when 
A = which coincide with the resonances of the potential (|59D found in Ref. The small 
open circles show their movement when A is gradually increased in steps of AA = 1 MeV. 
The digits displayed in Table III are those which do not change when the rotation angle 
changes and thus they represent the accuracy of our calculations. 

V. CONCLUSIONS 

The present work is a continuation of a series of papers 0,|^J31 which a practical 
method for quantum mechanical calculations is developed. The method is based on direct 
calculations of the Jost function and Jost solutions and is a combination of the variable- 
constant method with the complex coordinate rotation. Here we have extended it to include 
non-central potentials. 

The proposed method offers a unified way of treating bound, scattering, and resonance 
state problems. It is a powerful method which enables us to investigate the analytical prop- 
erties of the Jost function in the complex momentum plane. This opens up new possibilities 
in locating resonances. Within this method the Siegert wave functions can be properly nor- 
malised. Even sub-threshold resonances can be located which is a difficult task for many 
other methods. Moreover the formalism presented can be easily extended to complex values 
of the angular momentum £ and therefore Regge trajectories can also be located 0. 

The wave function can be obtained in a form which guarantees its correct asymptotic 
behaviour for all three types of physical problems. In all cases the same accuracy can be 
achieved which can be reliably controlled by simply changing the rotation angle. 

Despite the restriction (^, the potentials with Coulomb tails can be incorporated into 
the proposed method in a straightforward way as it was done in our previous publications 
To do this we need only to replace in all formulae the Riccati-Bessel functions by the 
corresponding Coulomb functions. 

The method can be extended further: Firstly, to treat the A^-body coupled hyperradial 
equations which differ from the two-body radial equations only by the possibility of having 
half-integer values of i; secondly, to investigate the behaviour of the coupled partial waves 
when the angular momenta is complex valued; thirdly, to treat non-analytical and singular 
potentials, and fourthly to treat coupled channel problems with different thresholds. The 
work on all these extentions is under way. 
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APPENDIX A: REGULAR BASIS 



By definition, ||$(A;,r)|| is a regular basis (fundamental matrix) if it has the following 
three properties fl]]: 

i) Each column of is a solution of Eq. (|Tl|) , 

ii) $[,,][,](A;,0) = 0, V[f ],[£], 

iii) all columns of ||<l'[^/][£]|| are linearly independent. 

Moreover, ||$|| must be a square matrix since Eq. (|TTp has as many independent regular 
column-solutions as the column dimension p4| . 



1. Behaviour at short distances 



The property ii) implies that the boundary condition (|T^) should be imposed at the 
point r = 0. Eq. (|ll]), however, is singular at the origin, and consequently the existence and 



uniqueness theorem [0] is not valid at this point. This theorem is valid for all r > 5 > 
for any arbitrary small 6. Hence, the matrix ||$|| can be uniquely defined by ||$(/c,(5)|| and 
c}j.||$(/c, 5) II, while within the interval [0, 6], \\^{k, r)\\ (having the above three properties) can 
be obtained explicitly as follows. After multiplying by and using (^, Eqs. ( pA]) decouple 
giving 



r^a^ - i{i + 1) (fc, r) ^ 0, re [0, 6] . 



£+1 



and 



e+1 



(Al) 

^ and 
and 



Each of these equations has two independent solutions behaving as 
of course the trivial (zero) solution. Therefore a regular column consists only of 
zero elements. Therefore the only way to construct the regular matrix ||$(A;, r)|| obeying 
Eq. (|A1|) with linearly independent columns, is to choose its diagonal as follows 



mk,r) 



( r^i+^ 

g ^^2+1 



\ 




V 







(A2) 



which obviously satisfies the above-mentioned three properties. Higher order corrections to 
Eq. ( [AID give, of course, non-zero off-diagonal elements. Each column of the matrix ( |A2| ) 
is a solution of Eq. (|A2|) and therefore has, separately, the short-range behaviour 



l(A;,r) \ 



I \ 




(A3) 
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Hence when r — all elements of the left hand side column of Eq. (^^) with [l\ ^ [^'\ 
vanish faster than the one with [£] = [£']. In other words, the off-diagonal elements of a 
column are infinitesimal as compared to the diagonal one. Therefore we can rewrite Eq. 
indicating the higher order corrections as follows 



/ r 



|$(/c,r)|| ~ 



+1^ 



,£2+1 



(A4) 



The inclusion of the higher order terms retains the linear independence of the columns. 
Indeed, 



det||$(fc,r)|| = n/^+i + |^n/'+^^ 



because the products involving off-diagonal elements always contain at least one of the 
factors o{r^^~^^), o{r^^~^^), ■ ■ ■, o{r^^~^'^). Hence 

det ||$(A;,r)|| — >TTr^'+i 

t 

and on the interval (0,5] we can always find at least one point tq where det ro)|| 7^ 0. 
This means that the columns of the matrix ( |A4| ) are linearly independent on the whole 
interval (0, 6]. 



The structure of the matrix ( |A4|) implies that the behaviour of the regular basis in the 
immediate vicinity of the point r = is such that 

Since we can choose the normalisation constant in each element of the matrix ||$(/c,r)|| 
independently, we may replace '^^ in the last equation by {krY +^ to get the boundary 
condition in the form 



r^o jt'[kr) 



mm 



(A5) 



To obtain the higher order terms of ( |A^ ) in an explicit form we replace ||$(fc,r)|| by the 
combination of the two new unknown matrix functions \\F^^\k^ r)||, as is given by Eq. (p!5[), 
which are subjected to the constraint (|l^) and obey Eqs. ([T7|) . Since by definition the 
matrix ||<I>(A;, r)|| is regular at r = 0, we have 



On the other hand 



hf\kr)F^^,^{k,r) -/i(-)(A:r)Fg,](A;,r) 

hg^\kr) ^ jejkr) + ini{kr) ^ _^ 

h[~\kr) ~ je{kr) - ine{kr) r-o 
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from which we conclude that in the immediate vicinity of the point r = the matrices 
r)|| and r)|| become identical, i. e., 

||FW(fc,r)|| — > \\A{k,r)\\, 

where the matrix ||A(A;, r)|| describes their common short-range asymptotic behaviour. 
Thus, we obtain 

$[^][^,](fc,r) — ^ ^ [h'i^\kr) + h^f ^kr)] A[^][e]{k,r) = je{kr)A[e^[e]{k,r) . (A6) 

Comparing the last two equations with Eq. ([A5D, we get the following boundary conditions 
for the matrices r)|| 

1™ . ; : = o\e]\e'] . ( A7) 

Note that unlike these conditions do not demand that all elements of the matrices 

be regular. Indeed, when £ > Eq. ( |X7| ) holds even if M^L(/c,r) ~ r^'y^^^'-<') with 



any e > 0. Hence, the left bottom corners of the matrices may, in principle, have 

diverging elements. As can be seen from the structure of the regular basis, Eq. (|15]), the 
functions r) are closely related to the Jost solutions and thus their singular be- 

haviour is not surprising. The r) itself remains always regular due to the presence 

of je{kr) in (|A6|) which compensates the diverging terms. 



Eq. ( |A7| ) gives us the boundary conditions in explicit form only for the diagonal elements 
of the matrices H-F'-^^H while for the off-diagonal ones it only implies that fI^^^^,-^ ~ o{ji/ je). 
To obtain them explicitly, let us take indefinite integrals in both sides of Eqs. ([171) . This 



gives 

1 f 

Fm']{k, r) = const ±— h^^\kr) J] (r)<l>[^»][£,](fc, r) dr . (A8) 

The integration constants are fixed for the diagonal and the right-top corner elements by 
the conditions ( |A7D . Indeed, due to (H) the integrals in (|A8|) for i < £' give functions 
j,-e^-2+e^e +1 6 > 0) which vanish at r = 0. Hence to fulfil ( [A7| ) we must set 
const=l for i = i' and const=0 for £ < £'. For the left-bottom elements we still have to 
choose the normalisation which is fixed by letting const = Thus finally 



p(±) 



1 r 

',^{k,r) = 6[i][fj] ± — / h'^P{kr)Y,W[e][en]{r)(^[e'][i']i^^^)dr, (A9) 



where the indefinite integrals should be understood as primitive functions since all arbitrary 
constants have been fixed already. 

Consider now Eqs. ([A9| ) on a small interval r G [0,5], where all functions under the 
integral can be replaced by their power-series expansions. While this is possible for hf^\kr) 
and iy[£][^/](r), for <I>[£][^/](A;, r) only the diagonal elements are known, 

<l>[^][<.](A;,r) ^ , re (0,5]. 
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As can be seen from Eq. ( |A4| ) the elements of the matrix ||$|| vanish (when r 0) with 
different speeds. In particular, in each row this speed increases from left to right. At 
the same time, within each column the element situated on the matrix diagonal has the 
lowest vanishing speed. This means that the leading term of the series expansion of a 
column is a column which is filled with zeros except the diagonal element. Looking at either 
the differential equations (0) or the integral equations (|A9|) , we see that in fact they are 
independent equations for each column. The series expansion of ||<l>|| should therefore be 
individually constructed for each column. Thus the leading term H^^^-*!! of such an expansion 
is 

= Ji{kr)6[e][i'] , re [0,6], (AlO) 
which, according to (|A6|) , implies that 

Fl^^^\k,r) = 5iem, re [0,6]. (All) 



Substituting ( |A10|) into the indefinite integral ( |A9|) and using the series expansions of hi, 



ji, and we obtain the first iteration || for all elements of the matrices 

and thus we can get H'J*'-^^!! which includes the next terms of the expansion of the regular 
solution. Using this iterative procedure we can find as many terms of the expansion as 
needed by the following recurrence formulae 



F{§!-^'\k,r) = ± ^ / hi^\kr)Y.W[e]inir)'^Pw']ik,r)dr, (A12) 



^\fak,r) = ;- hr>{kr)F,^>\-\k,r) + hr\kr)Fl,^^^^ . (A13) 



me'] 



2. Behaviour at large distances 

To analyse the long-range behaviour of the matrix-functions ||-F*^^^(/i;, r)|| we rewrite the 
system ([T7| ) as follows 

a,F(^)](fc,r) = ± ^h^r\kr) ^W^[,][,.](r)$[,„][,,](A;,r) (A14) 

It is clear that if the limits lim ||F*^^)(A;, r)|| exist, then the corresponding derivatives 

5r||-F*^^''(A;, r)|| are zero at large distances. The converse is valid, however, only if the deriva- 
tives vanish faster than r~^. Indeed, when the derivative of a function (f{r) behaves as 

the asymptotic behaviour of the function itself can be written as 

Inr e = 



(p{r) ~ J r ^^^^^dr = const + < 1 



e^O 
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which obviously has a finite hmit only if e > 0. 



Therefore, to establish the domains of the fc-plane where the limits lim llF'^^-'fA;, r)|| 

r— »oo 

exist, it is sufficient to find out at which areas of this plane the right-hand sides of Eqs. 
( |A14| ) vanish at large distances faster than the Coulomb potential. For this we, first of 
all, need to know how the basic solutions (/c, r) behave in different domains of the 

fc-plane when r —>■ oo. To this end we consider Eqs. (|TT|) at large r. If the potential is of 
short-range (decaying exponentially or faster) these equations are reduced to the un-coupled 
Riccati-Bessel equations, viz.. 



0. 



(A15) 



Considering them as one matrix equation, then there are two linearly independent solutions 
which can be chosen as 



\\h^^\kr)\\ = 



I hf\kr) 








(A16) 



V ■■■hfJ{kT)j 

Any other solution is a linear combination of them. Therefore 

mk,r)\\^\\h^^\kr)\\.\\C4 + \\h^-\kr)\\^\C, 
and according to (|2^) 

||$(fc,r)|| — > exp(?fcr)||C'i|| + exp(--iA;r)||C'2|| . 



(A17) 



(A18) 



If Im kr = Q then both terms of the last asymptotic relation are oscillating and therefore 
the regular solution is bounded, i.e.. 



{k,r] 



< const, V 



when r oo 



(A19) 



For Imfcr 7^ 0, the regular solution diverges due either to the first or the second term of 
( |A18|) . In this case we can keep only the diverging term of it, viz.. 



||$(A;, r)|| ^~^exp(-iA;r)||C2| 
mk,r)\\^^expi+tkr)\\C,\ 



Im A;r > , 
Im fcr < . 



(A20) 
(A21) 



There are, however, special points ko on the A;-plane, called spectral points, where by defi- 
nition the regular solution contains only the first term of ( |A17| ), i.e.. 



<i>{ko,r)\\^ exp(+iA;or)||Ci|| . 



(A22) 



Spectral points situated on the positive imaginary axis correspond to bound states, and 
those situated under the real axis to resonance states. 
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Though we have derived Eqs. ( |A19| , [A20| , [A21| ) describing the asymptotic behaviour of 



the regular basis, under the assumption that the potential is of the short range, they should 
be valid also for potentials decaying faster than ~ 1/r. The general proof, however, requires 
more sophisticated mathematical methods which are beyond the scope of this article. 



Substituting Eqs. ([M9| ), (|X20D , and (1X22|) into (^T|) and using the asymp 



totic form of the Riccati-Hankel functions, Eq. (^), together with the constraint on the 
long-range behaviour of the potential, Eq. (^, we find that in general the derivatives of 
||F'^+)(/c, r)|| and r) || vanish rapidly enough in the different domains of the fc-plane. 

The only points where both limits lim r)|| exist simultaneously are those with 

Im fcr = and the spectral points with Im kr > 0. 

For real r then lim IIF'^+HA;, r)|| exists for ImA; < and at the spectral points which 

r— +00 

correspond to bound states, while lim ||F'^~)(A;, r)|| exists for ImA; > 0. In this case the 

1 — ►oo 

dividing line that separates the two domains of the fc-plane coincides with the real axis. 
This line can be turned downwards, to expose the resonance spectral points, by rotating r 
as given by Eq. (P7|). Indeed, if (p is the polar angle parametrising a point on the fc-plane 

k = \k\e'^, 

then by choosing large enough > we can make Imkr, 

ImA;r = Im (|A;|xe^('^+^)) = |A;|x sin(^ + 0) , (A23) 



positive even when is negative, —6 < (p < tt — 6, i.e., when the point k is in the fourth quad- 
rant. From the last equation is clear that when 6^0 the dividing line is 
(see Fig. p. 



It is worthwhile to mention that the border separating the two domains of the complex Ac- 
plane is a line only in the case of long-range potentials obeying the condition (^. If however 
the potential decays at large r exponentially, then lim ||F*^+)(A;, r)|| exists also within a band 

above the dividing line while lim ||F(~)(A;, r) || in the symmetric band below this line. The 
faster the potential decays the wider this band is. For example, if the potential decays as 
~ exp(— /ir) then the right hand side of the second equation of the system ( |A14| ) behaves, 
below the dividing line (ImA;r < 0), as 

drF^-\k,r) ~ e^'^'e-'^'e^'^" = e^(2Refcr-/.imr) (_2lmA;r - /iRer) . 
The last exponential function decays at large distances if ImA;r > — ^Rer, i.e., when 

|A;|sin(^ + 0) > -^cos^ . (A24) 

If 6 = 0, this condition reads ImA; > — yu/2 which enables us to locate also virtual states 
(spectral points on the negative imaginary axis) situated not far from the origin. 
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TABLES 



TABLE I. Bound states generated by the RSC and Moscow potentials 





our method 


Refs. 




potential 


ko (fm-i) 


Eo (MeV) 




Eo (MeV) 


D% 


RSC 


iO. 2316110 


-2.22460 


6.470 


-2.22460 


6.470 


Moscow 


i3.5571773 


-524.741 


14.36 


-524.8 




Moscow 


i0.2316000 


-2.22439 


6.588 


-2.2246 


6.778 



TABLE II. Jost matrix for the RSC potential at A: = 0.5 exp(— z0.37r) fm , calculated with 
Xmax = 50 fm and three different values of the rotation angle 6. 



e 


\\F^-\k,Xma.e'')\\ 





( 22637458292890 + il4763693403731 -22818484 + ^84501579 \ 
y 8536691059869 + i6193286881927 -10532586 + ^32492267 / 


0.357r 




^ -325291222087 + i639117048997 -2294097 - i362305 ^ 
^ 145085673415 - ^1488370181212 4729536 - il042674 ) 




0.407r 




' -325288820471 + i639116901481 -2294097 - i362305 ^ 
^ 145081315229 - ^1488368045635 4729536 - il042675 ) 
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TABLE III. Spectral points of the model potential (see text) with A = 15 MeV. 



Kepo (fm ^) 


Impo (f™ ^) 


Eo (MeV) 


r (MeV) 





4.5581531714 


-10.3883801672 








4.023079863 


-8.09258579 








3.471206989 


-6.02463898 








2.899849780 


-4.20456437 








2.305427868 


-2.65749883 








1.681898798 


-1.41439178 








1.01619035 


-0.5163214 








0.254097 


-0.0322826 





Q /1/|fi«nSQ9 
0.44DDUoyZ 


— U.0oUii4oy 


c 7QQn/i i^Qn 


q /1 1 Q/in 
O.D04iy4U 


4.1388078309 


-0.1467148600 


8.5541025053 


1.214449223 


4.46521096 


-0.686071761 


9.73370724 


6.12691030 


-1.71-1321 


-1.3323(35 


10.36071 


12. (3 1231 


4.96356 


-1.99719 


10.3241 


19.8263 


5.1410 


-2.6634 


9.6681 


27.385 
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FIGURES 
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FIG. 2. Spectral points (filled circles) of the model potential with A = 15MeV. The large 
open circles represent three of the resonances generated by the central potential Vc and the small 
circles show the movement of the resonances when A decreases from 15MeV to zero in steps of 
AA = 1 MeV. The dividing line corresponds to ^ = 0.27r. 
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